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Abstract. In this work the finite /3-cffccts of an electron branch of the geodesic 
acoustic mode (el-GAM) driven by electron temperature gradient (ETG) modes is 
presented. The work is based on a fluid description of the ETG mode retaining 
non-adiabatic ions and the dispersion relation for el-GAMs driven nonlinearly by 
ETG modes is derived. A new saturation mechanism for ETG turbulence through 
the interaction with el-GAMs is found, resulting in a significantly enhanced ETG 
turbulence saturation level compared to the mixing length estimate. 
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1. Introduction 

Experimental investigations has elucidated on the complex dynamics of the low to high 
(L-H) plasma confinement mode transition. Evidence of interactions between the the 
turbulence driven E x B zonal flow oscillation or Geodesic Acoustic Mode (GAM) 
PQ|2i|3j|Il[7J[6j[EJ|9], ^ ur ^ u ^ ence an< ^ the mean equilibrium flows during this transition 
was found. Furthermore, periodic modulation of flow and turbulence level with the 
characteristic limit cycle oscillation at the GAM frequency was present [2]. The GAMs 
are weakly damped by Landau resonances and moreover this damping effect is weaker 
at the edge suggesting that GAMs are more prominent in the region where transport 
barriers are expected [3]. 

During recent years investigations on coherent structures such as vortices, streamers 
and zonal flows (m = n = 0, where m and n are the poloidal and toroidal modenumbers 
respectively) revealed that they play a critical role in determining the overall transport in 
magnetically confined plasmas |4J . Zonal flows impede transport by shear decorrelation 
[U [5] , whereas the GAM is the oscillatory counterpart of the zonal flow (m = n = in 
the potential perturbation, m = 1, n = in the perturbations in density, temperature 
and parallel velocity) and thus a weaker effect on turbulence is expected E]. 

The electron-temperature-gradient (ETG) mode driven by a combination of 
electron temperature gradients and field line curvature effects is a likely candidate 
for driving electron heat transport [TOj [TTJ [T2J [T3J [H] . The ETG driven electron 
heat transport is determined by short scale fluctuations that do not influence ion heat 
transport and is largely unaffected by the large scale flows stabilizing ion-temperature- 
gradient (ITG) modes. The generation of large scale modes such as zonal flows and 
GAMs is here realized through the Wave Kinetic Equation (WKE) analysis that is 
based on the coupling of the micro-scale turbulence with the GAM through the WKE 
under the assumptions that there is a large separation of scales in space and time. 

Although the electron GAM (el- GAM) was previously explored in Refs |2T] and 
|22j, vital physics were however neglected, in the present work the finite /3-effects are 
elaborated on as well as numerical quantifications of the frequency and growth rate are 
given. The finite /3-effects are added in an analogous way compared to the recent work on 
zonal flows in Ref . [23| 121] . In particular the Maxwell stress is included in the generation 
of the el-GAM. The frequency of the el-Gam is higher compared to the ion GAM by 
the square root of the ion-to-electron mass ratio (fl g (electron)/Q q (ion) ^mi/m e 
where Q q (electron) and Q q (ion) are the real frequencies of the electron and ion GAMs, 
respectively). 

2. The linear Electron Temperature Gradient Mode 

In this section we will describe the preliminaries of the electron-temperature-gradient 
(ETG) mode which we consider under the following restrictions on real frequency and 
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wave length: fij < u ~ u± « f2 e , k±Ci > u > kuc e . Here fL are the respective 



cyclotron frequencies, pj the Larmor radii and Cj = ^Tj/mj the thermal velocities. 
The diamagnetic frequency is u± ~ kgp e c e /L n , k± and k\\ are the perpendicular and the 
parallel wave vectors. The ETG model consists of a combination of an ion and electron 
fluid dynamics coupled through the quasineutrality including finite /3-effects [131 EH] ■ 

2.1. Ion and impurity dynamics 

In this section, we will start by describing the ion fluid dynamics in the ETG mode 
description. In the limit u > k\\c e the ions are stationary along the mean magnetic 
field B (where B = B e\\) whereas in the limit k±Ci » u, k±pi » 1 the ions 
are unmagnetized. In this paper we will use the non-adabatic responses in the limits 
uj < k±ci < k±Ci, ci = an d assume that fij < u < Vt e are fulfilled for the ions and 
impurities. In the ETG mode description we can utilize the ion and impurity continuity 
and momentum equations of the form, 
dn ' 

— l + n j V-v j =0, and (1) 

at 

dv ' 

m jnj — + e n i V0 + T i Vn i = O, (2) 

where j = i for ions and j = I for impurities. Now, we derive the non-adiabatic ion 
response with r« = T e /Tj and impurity response with with 77 = T e /Tj, respectively. We 
have thus, 

= - ( i-J/Wff ) ~ (3) 

Here Tj and nj are the mean temperature and density of species (j = e,i,I), where 
hi = Sn/rii, hj = 5nj/ni and = e<p/T e are the normalized ion density, impurity 
density and potential fluctuations. Next we present the electron dynamics and the 
linear dispersion relation. 

2.2. The electron model 

The electron dynamics for the toroidal ETG mode are governed by the continuity, 
parallel momentum and energy equations adapted from the Braginskii's fluid equations. 
The electron equations are analogous to the ion fluid equations used for the toroidal 
ITG mode, 

+ V • (n e v E + n e v* e ) + V • {n e v pe + n e v m ) + V • (n e v\\ e ) = (4) 
3 30T 

-n e —^- + n e T e V-v e + V-q e =0. (5) 

Here we used the definitions q e = —(5p e /2m e fl e )e\\ x VT e as the diamagnetic heat 
flux, ve is the E x B drift, v+ e is the electron diamagnetic drift velocity, vp e is the 
polarization drift velocity, v n is the stress tensor drift velocity, and the derivative is 
defined as d/dt = d/dt + p e c e e x V0 • V. A relation between the parallel current density 



High Frequency Geodesic Acoustic Modes in Electron Scale Turbulence 



4 



and the parallel component of the vector potential (Jn) can be found using Ampere's 
law, 

Vii|, = -^J,,, (6) 

Taking into account the diamagnetic cancellations in the continuity and energy 
equations, the Eqs. OH [5] and [6]) can be simplified and written in normalized form 

as, 

- w - vi Jr^ - ( x + (1 + %)vl ) - V|1 vii|1 + 

( Id d \ ~ 

-(/3 e /2)[i||,Vji||] + [0,V 2 0], 

K/3 e /2 - Vi)^ + (1 + r /e )(/3 e /2)V^ A,, + V N (0 - n e - f e ) 

72 



-(/3 e /2)[0-n e ,A,|] + (/3 e /2)[T e ,A||] + [0, Vi^,], 
-T e + -e n costf-- + sm^ -_T e + fo e - -)--* - -- 



Note that similar equations have been used previously estimating the zonal flow 
generation in ETG turbulence and have been shown to give good agreement with linear 
gyrokinetic calculations [131 EH]- The variables are normalized according to 

(0,n,T e ) = (L n /p e )(e5(fr/T eo ,5n e /n ,5T e /T e0 ), (10) 

i,l = {2c e L n /p e c Pe )eA ll /T e0 , (11) 



P e = 8nnT e /B 2 . 

Using the Poisson equation in combination with (|3]) we then find 



n, 



+ 



+ k 2 X 2 



De 



(12) 



(13) 



\l-C0 2 /k 2 ± C 2 l-U 2 /(k 2 ± C 2 , 

First we will consider the linear dynamical equations (0 El and |9]) and utilizing Eq. ([13 
as in Ref. [H] and we find a semi-local dispersion relation as follows, 



(jj- 



A e + ^(1 + A e ) ] + (1 - e„(l + A e )) 



k ±Pe (W - (1 + 



+ 
+ 



) ( (r? e - ^)w* + ^wA f 



2 7,2 7,2 2 



:i+a £ 



(14) 



In the following we will use the notation A e 



Ti(ni/n e )/{l - u 2 /k\c 2 ) + 



r,(Z e// n,/n e )/(l - co 2 /k 2 ± c 2 ) + fc 2 A 2 



Note that in the limit T,- 
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k±\oe < k±p e < 1 and in the absence of impurity ions, A e « 1 and the ions follow the 
Boltzmann relation in the standard ETG mode dynamics. Here Xoe = \Jt c / (4:iin e e 2 ) is 
the Debye length, the Debye shielding effect is important for Xoe/Pe > 1- The dispersion 
relation Eq. (TH]) is analogous to the toroidal ion-temperature-gradient mode dispersion 
relation except that the ion quantities are exchanged to their electron counterparts. 



3. Modeling Electron Geodesic Acoustic modes 

The Geodesic Acoustic Modes are the m = n = 0, k r ^ perturbation of the potential 
field and the n = 0, m = 1, k r ^ perturbation in the density, temperatures and 
the magnetic field perturbations. [U 0] The el-GAM (q, Q q ) induced by ETG modes 
(k, u) is considered under the conditions when the ETG mode real frequency satisfies 
Q e > u > Qi at the scale k±p e < 1 and the real frequency of the GAM fulfills Q q ~ c e /R 
at the scale q r < k r . 



3.1. Linear Electron Geodesic Acoustic Modes 

We start by deriving the linear electron GAM dispersion relation following the outline 
in the previous paper Ref . [21] , by writing the m = 1 equations for the density, parallel 
component of the vector potential, temperature and the m = of the electrostatic 
potential 

" + e n sin6^ - V„Vii« = 0, (15) 

(/3 e /2-Vi)|i«-V||(n2 + T«) =0, (16) 

Hf (1) - -— - o fin 



-V^4 0) -^sin^(n2 + f«) =0. (U 



We will derive the linear GAM frequency as follows by using Eq. ( f!6|) and by eliminating 
the density m — 1 component using a time derivative of Eq. f llSp .Finally by utilizing 
Eq. (US]) we find, 



Pe^Vi0(°) + e n v*{sin9^ L^sin^ + V„i- )) = 0. (19) 



Here (■ • •) is the average over the poloidal angle 6. In the simplest case this leads to the 
dispersion relation, 

°Hi( 2+ ?TT«k!))- (20) 

Note that the linear electron GAM is purely oscillating analogously to its ion counterpart 
c.f. Refs. [6]. Here it is of interest to note that it is very similar to the result found in 
Ref. [221. 
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3.2. The Non-linearly Driven Geodesic Acoustic modes 

We will now study the system including the non-linear terms and derive the electron 
GAM growth rate. The non-linear extension to the evolution equations presented 
previously in Eqs © - ([9} are 



()f Vl-0 - (1 + (1 + Vem) V e <P - V||V1A|| + 

( Id d\ ~ 

€n { C ° S 9 rd9 + 9 dr~ ) ^ ~ ™ e ~ Te) 

+ [0,V 2 0]-(/3 e /2)[i||,V 2 i||], 

f(/3 e /2 - Vi)^ + (1 + r] e )(/3 e /2)V e ) A\\ + V N (0 - n e - f e ) 



-T e + -e n loose-- + an*- -_T e + fo e - -)--* - --n e 



In order to find the relevant equations for the electron GAM dynamics we consider the 
m = 1 component of Eqs (12TT) - ( 1231) . 



-^ + e " Sm ^ 0) - V '' Vii ^ = 

([0 fe ,V 2 fc ])«-(/3 e /2)([I|| fc ,V 2 i|| fc ])« =0, (24) 
(&/2 - V 2 ± )|i« - V||(n2 + T«) = 

([0 fc ,Vli||,]) (1) =0, (25) 

^_T (1) -~-fi (1) 
W eG 3di HeG 

-{[kXk]) {V) =N[ l \ (26) 

where superscript (1) over the fluctuating quantities denotes the m = 1 poloidal mode 
number and (■ ■ •) is the average over the fast time and spatial scale of the ETG 
turbulence and that non-linear terms associated with parallel dynamics are small since 
■%<<!. We now study the m = potential perturbations, 



V72 ^ 7(0) • Q d 

([fa, V 2 fc ]) (O) - W2)([A\\ k , V 2 i|| fc ])W = iV 2 (0) . (27) 

In order to evaluate the Maxwell stress part in Eq. fT271) , we will approximate the parallel 
part of the electromagnetic vector potential with the electrostatic potential through a 
linear relation. The relation Ak\\ = Aofa is found by using the Eqs. ([8]), ()9]) and the 
non-adiabatic response Eq. (fl3l) giving an approximation of the total stress of the form, 

Arf = (l-|fi a | 2 )([4,V 2 4])(°). (28) 
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The fl a factor is found by using Eq. (jSJ) 



I a\ 2 



(29) 



(P e /2 + kl)u-(l + Ve )(3 e k e /2 
where (po is determined by the temperature equation 
~ (r ?e -2/3)-2/3a;A e I ~ 

Tek = -^75 ; 4>k = <fo<Pk, (30) 

u + 5/3e n gk y 

and A e is determined by the non-adiabatic response condition. Note that fl a vanishes 
at /3 e — 0. In the above we defined non-linear term on the RHS in Eqs. ( 1241) - (1271) 
as = J dtp\c e z x V0 • VV^_0. This can be written T e = + A"^, where the 
m = 1 component is determined by an integral of the convective non-linear term as 
Nl = — J dtp s c s z x V^ - 1 • VT e W. This leads to a relation between the m = 1 component 
of the density and temperature fluctuations modified by a non-linear term. Here, the 
non-linear terms can be written in the form, 

^ (1) =E^V r |0 fc | 2 , (31) 
iV 2 (0) = {I - \n a \ 2 )qlY.Kke\k\ 2 ■ (32) 

k 

We continue by considering the Equations (124"]) and (1271) for the m = 1 component 
and m = component, respectively, 

fL^_ e ^fUa> (33) 

at en or 
^Vim + e n (s^l^) = N^. (34) 

We keep the A"} 1 ^ non-linear term in order to quantify the effects of the convective non- 
linearity. Similar to the operations performed to find the linear electron GAM frequency 
we eliminate the m = 1 component of the electron density by taking a time derivative 
of Eq. (134]) this yields, 

^Vi4 01 4 + < + |<)> = (35) 

Note that the el-GAM wave equation will be modified by the effects of the parallel 
current density (J\\) and the m = 1 non-linear terms in the general case, however we 
see by inspection that on average the term A"} 1 ^ does not contribute whereas the 
non-linearity may drive the GAM unstable. We will use the wave kinetic equation 
HI H51 HSl HH El HS1 HSl ED] to describe the background short scale ETG turbulence 
for (fl q ,q) < (u),k), where the action density N k = E k /\u r \ ~ eol^fcl 2 / ^- Here e |</>fc| 2 , 
is the total energy in the ETG mode with mode number k where eo = t + k\ + -rji- 
The electrostatic potential is represented as a sum of fluctuating and mean quantities 
<f)(X,x,T,t) = $(X,T) + <p(x,t) where <3>(X,T) is the mean flow potential. The 
coordinates (x,T^, (x,t) are the spatial and time coordinates for the mean flows and 
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small scale fluctuations, respectively. The wave kinetic equation can be written as 



= 7*JV*(r, «)-AwJV t (r,«) 2 . (36) 



We will solve Equation ( |36|) by assuming a small perturbation (5N k ) driven by a slow 
variation for the GAM compared to the mean (N k o) such that N k = N k o + 5N k . The 
relevant non-linear terms can be approximated in the following form, 

v 2 Jk\) ~ (i - \n a WrJ2 k rke— SN k (q,n q ). (37) 

k e o 

For all GAMs we have q r > qe, with the following relation between d~N k and dN k o/dk r , 

r 71 r • 2, ,0 p ^Ofc , kgq r f^N ok 

6N k = -iq r k e cj) G R—— H (38) 

C fc r T iy J (Tie - 7}^) 

where we have used 5u; 9 = ■ ~ i(kgq r — k r qe)(j)Q in the wave kinetic equation and 

the definition R = „ — - — . Using the results from the wave- kinetic treatment we 

can compute the non-linear contributions to be of the form, 



-i(l-\n^)q^k r kV-^R 



e ok r 



+ (1-|^| 2 )^E^- 7 — u /2 "2- (39) 

In order to find the non-linear growth rate of the electron GAM we need to find relations 
between the variables n^l, T^q and <j>Q, 

«3 = - ^ 3i ",v' it (40) 

Q2 _ 5 ^ 

Using Eq. (140]) and the Fourier representation of Eq. (J35l) resulting in, 

gr 9 3/3 e /2 + g r 2 9 6 n9r 9 



;i - i^i 2 ) k - f^S-^ (4i) 



v " 9 3 /3 e /2 + g r 2 ; ^ ^ e " «9fc r 

Eq. (HTj) is the sought dispersion relation for the electron GAM. In the electrostatic 
limit we find a perturbative solution of the form, 



7<? _ lqfpikope 1 



7 -^n 

0fc 

Pe 



(42) 



Here the main contribution to the non-linear generation of el-GAMs originates from 
the Reynolds stress term in competition with the Maxwell stress term. Note that the 
Maxwell stress term reduces the drive of the GAM analogously to the zonal flow situation 
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Figure 1. (color online) The el-GAM growth (blue, solid line), real frequency (blue, 
dashed line) and the linear solution (red, dash-dotted line) as a function of r\ e is 
shown for the parameters e n = 0.909, /3 — 0.01 k r — 0.6 = kg and q = 1 in the 
strong ballooning limit g(9) = 1 and k\(9) = k% + kg with ETG saturation level 
\4>k\ ~ {lk/u*){l/k x L n ), 




Figure 2. (color online) The el-GAM growth (blue, solid line), real frequency (blue, 
dashed line) and the linear solution (red, dash-dotted line) as a function of g e is 
shown for the parameters g e — 4.0, e n = 0.909, f3 — 0.01 k r = 0.6 = kg in the 
strong ballooning limit g(0) = 1 and k\(9) = k% + kg with ETG saturation level 

\4>k\ ~ (jh/u-Jil/kxLn). 



[231 and that this results differ to the result found in EH where Q n = 0. The non 



linearly driven electron GAM is unstable with a growth rate depending on the saturation 

~ 2 

level <pk of the ETG mode turbulence. 

We will now present the numerically determined frequencies and growth rates found 
by solving Eq. ( HIT) in few limits. We have chosen common parameters relevant for an 
ETG mode situation with r\ e = 4.0, e n = 0.909, (3 = 0.01 k r = 0.6 = k g and q = 1 in the 
strong ballooning limit g(6) = 1 and k\{9) = + kg. Furthermore the saturation level 
of the ETG mode is represented by the mode coupling estimate « ( , yk/w*)(l/k x L n )- 



High Frequency Geodesic Acoustic Modes in Electron Scale Turbulence 



10 



In Figure 1, we display the el-GAM Q q growth rate and real frequency (in blue) as 
a function of r] e in comparison with the solution to the linear el-GAM dispersion relation 
Eq. ()20lin red). The el-GAM growth rate is increasing with increasing t] e analogously 
to the linear ETG growth rate due to the mode coupling saturation level. 

Figure 2 shows the el-GAM Q q growth rate and real frequency (in blue) as a 
function of q in comparison with the solution to the linear el-GAM dispersion relation 
Eq. f)20l in red). It is found that q has a stabilizing effect on the el-GAM however the 
quantitative results are strongly dependent on the parameters while the finite ^-effects 
are moderate. The effect of the Maxwell stress on the growth rate is rather small on the 
order of (3 e ~ 0.01 for the given parameters however for increasing amount of impurity 
ions (increasing Z e ff) it becomes increasingly important. Note that, in the dispersion 
relation a set-off non-linear drive is present below which the el-GAM is stable. 



4. Saturation mechanism 



In this section we will estimate a new saturation level for the ETG turbulent electrostatic 
potential (<f) k ) by balancing the Landau damping in competition with the non-linear 
growth rate of the GAM in a constant background of ETG mode turbulence, according 
to the well known predator-prey models used, c.f. Eq. (4) in Ref. [7], 

dN k 



dt 

dU G 
dt 



lk N k - AujNI 2 - ^U G N k 
IqUc ~ 1lU g - v*U G . 



(43) 
(44) 



Here we have represented the ETG mode turbulence as N k 



| 2 § and Uq 



sin8) with the following parameters: 7 is the ETG mode growth rate, 71 
is the coupling between the ETG mode and the GAM. The Landau damping rate 



(~ T = £s. 

the neoclassical damping in stationary state 



is assumed to be balanced by GAM growth rate Eq. (|42|) modified by 

— > and — > 0. In steady state find 
1l + v% 



8N 
dt 



the saturation level for the ETG turbulent intensity as (p/ q 



e<f) k L r 



Pe 



2L n 
~q~R 



2^2 ) V^Ve\ -\\ Z 




(45) 



Here, the saturation level is modified by the neoclassical damping v* = 
— factor arises due to the spatial extension of the GAM and we obtain, 

q r 



v e ^- and the 

Veth 



e(f>k L r 



Pe 



30 - 40. 



Note that the result found using a mixing length estimate with 
smaller. Here in this estimation we have used L n = 0.05, q = 
l/q r ~ (pgLr) 1//3 , kgp e = 0.3 where kg/q r ~ 4 and r/ e ~ 1. 



e</> L n 

3.0, 



(46) 

~ 1 is significantly 
R = 4, e„ = 0.025, 
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5. Conclusions 

In this work the electromagnetic effects on the electron Geodesic Acoustic Mode 
(el-GAM) are investigated. The work extends a previous study (Ref. [2TJ) by self- 
consistently including linear as well as nonlinear /3 effects in the derivation. The linear 
dispersion relation of the el-GAM is purely oscillatory with a frequency Q q ~ ^ whereas 
the GAM growth rate, is estimated by a non-linear treatment based on the wave-kinetic 
approach. The el-GAM growth rate is driven by a competition between the Reynolds 
stress and the Maxwell stress. The non-linear dispersion relation is solved numerically 
where it is found that the magnetic safety factor q has a stabilizing effect on the el-GAM 
however the quantitative results are strongly dependent on the other physical parameters 
while the finite /3-effects are moderate. The effect of the Maxwell stress on the growth 
rate is rather small on the order of /3 e ~ 0.01 for the given parameters however for 
increasing amount of impurity ions (increasing Z e ff) it becomes increasingly important. 
Note that, in the dispersion relation a set-off non-linear drive is present below which 
the el-GAM is stable. 

To estimate the ETG mode fluctuation level and GAM growth, a predator-prey 
model was used to describe the coupling between the GAMs and small scale ETG 
turbulence. The stationary point of the coupled system implies that the ETG turbulent 
saturation level <pk can be drastically enhanced by a new saturation mechanism, 
stemming from a balance between the Landau damping and the GAM growth rate. 
This may result in highly elevated particle and electron heat transport, relevant for the 
edge pedestal region of H-mode plasmas. 
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